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^q' The ordering of the three-dimensional Heisenberg spin glass with the weak random anisotropy 

in magnetic fields is studied by extensive equilibrium Monte Carlo simulations. Both the spin and 
the chirality are monitored. We find strong numerical evidence that a replica symmetry breaking 
transition occurs in the chiral sector, which accompanies the simultaneous spin-glass order. Despite 
the similarity in the global symmetry, the ordering behavior of the weakly anisotropic Heisenberg 
spin glass differs significantly from that of the strongly anisotropic Ising spin glass. The obtained 
phase diagram in the temperature - magnetic field plane is similar to the experimental phase diagram. 
Our results highlight the importance of the chirality in the spin-glass ordering of the Heisenberg-like 
spin glass, and support the spin-chirality decoupling-recoupling scenario of spin-glass transitions. 



I. INTRODUCTION 



Spin glasses (SGs) are random magnets where ferromagnetic and anti-ferromagnetic exchange interactions co- 
exist and competed Experimentally, it is now well established that SG magnets in zero field exhibit a thermody- 
1 . ■ namic second-order phase transition at a nonzero temperature into the thermodynamic SG phase. By contrast, 
whether SG magnets exhibit a thermodynamic phase transition in applied magnetic fields has been a long-standing, 
q yet unsolved issue. This issue is closely related to the fundamental question of whether the SG ordered state in 
' zero field accompanies an ergodicity breaking not directly related to the global symmetry of the Hamiltonian, i.e., 
the replica symmetry breaking (RSB). 

The experimental evidence of an in-field transition of SG remains to be obscure. For the strongly anisotropic 
Ising-like SG, Feo.sMno.sTiOa, the non-existence of an in-field SG transition was reported in Ref.«. Meanwhile, 
■ many of real SG materials are more or less Heisenberg-like rather than Ising-like in the sense that the magnetic 
C^) ' anisotropy is considerably weaker than the isotropic exchange interaction* 1 Recent experiments on such weakly 
anisotropic Heisenberg-like SGs suggested the occurrence of an in-field SG transition; 3 ^ in apparent contrast to 
j^j ' Refi^. Setting aside the question of the strict nature of the apparent SG "transition" observed experimentally 
in applied fields, it has been known that the "transition line" between the paramagnetic and the SG phases is 
similar to the one obtained for the mean- field Sherrington-Kirkpatrick (SK) model: 5 Namely, in the weak field 
regime, the in- field transition temperature, T g {H), is rapidly suppressed with increasing the field intensity H, as 
H oc |1 - T g (H)/T g (Q)\ 3 / 2 (dc Almcida-Thouless (AT) line 6 ), while in the high field regime, T g (H) stays rather 
robust against H, behaving as H oc |1 - T s (7J)/T g (0)| 1/2 (Gabay-Toulouse (GT) line 7 .). However, the reason why 
the mean-field results have given such a good description of the phase boundary, including the values of the critical 
""2 ■ exponents describing the phase boundary, has remained to be a mystery. 

On the theoretical side, most of numerical studies on the finite-ranged SG models in three dimensions (3D) 
have focused on the properties of the Ising SGi Since no global symmetry exists in the Ising SG under magnetic 
fields, an in-field transition, if any, should be a pure RSB transition. Unfortunately, numerical simulations on the 
Ising SG have been unable to give a definitive answer concerning the existence of a thermodynamic SG transition 
in magnetic fields, 8 . 9 . 10 ! 11 . 12 ! 13 . 14 ! 15 ! 16 ! 17 ! 18 
$_i ' For the isotropic 3D Heisenberg SG in zero field, it has been believed for years that the SG transition occurs 
?3 i only at zero temperature, i.e., Tsq = Q^S&SL^S^ Since applied magnetic fields make the SG transition even more 
unlikely, one has expected no phase transition to occur in applied fields, and hence, until quite recently, no extensive 
numerical simulation has ever been performed for the 3D Heisenberg SG in magnetic fields. Meanwhile, recent 
studies have revealed that the Heisenberg SG possesses an important physical ingredient absent in the Ising SG, i.e., 
the c/uroHfa/i 24 ! 25 ! 26 ! 27 ' 28 ! 29 ! 30 ! 31 ! 32 In the chirality scenario of Refi 24 i 25 , in particular, the chirality is claimed to be a 
hidden order parameter of the SG transition of real Heisenberg-like SG magnets: In the fully isotropic Heisenberg 
SG, the spin and the chirality, though they are coupled at short length scales, are eventually decoupled at long 
length scales, and the system exhibits a chiral-glass transition at a finite temperature without accompanying the 
standard SG order. The chiral-glass transition corresponds to the spontaneous breaking of the spin-reflection 
symmetry with preserving the 5*0(3) spin-proper-rotation symmetry. In the more realistic weakly anisotropic 
system, the Heisenberg spin, decoupled from the chirality in the isotropic system, is "recoupled" to the chirality 
at long length scales via the random magnetic anisotropy. The SG order of the weakly anisotropic Heisenberg SG 
is then dictated at long length scales by the chirality ordering of the isotropic system. Some numerical support of 
such a spin-chirality decoupling-recoupling scenario was already reported in zero field ^^S^L^ 



By contrast, some other groups claimed that the chiral-glass transition of the 3D isotropic Heisenberg SG alread$ 
accompanied the standard SG order, which means that the standard SG order occurs at a finite temperature 
simultaneously with the chiralityii&iMiSiSi Note that this is in contrast to the earlier belief in the community 
that the SG transition occurs only at T = in the 3D Heisenberg SG. Ref. 29 maintains, however, that the SG 
order occurs at a temperature lower than the chiral-glass transition temperature, i.e., Tqg > ^SG > 0, and the 
controversy remains. 

Recently, the present authors performed the first extensive MC simulation of the 3D Heisenberg SG in magnetic 
fields, and have observed that the chiral-glass transition, essentially of the same character as in the zero-field one, 
occurs at a finite temperature even in magnetic fields^2*2i The chiral-glass transition line in the temperature- 
magnetic field phase diagram turned out to have a striking resemblance to the GT-line observed experimentally, 
although the nature of the transition is entirely different from the mean-field GT-line. Note that the fully isotropic 
Heisenberg SG in fields possesses the global Z 2 x 5*0(2) symmetry, the chiral Z 2 referring to the global spin- 
reflection with respect to the plane containing the magnetic-field axis, and the SO(2) referring to the global 
spin-rotation around the magnetic-field axis. 

In the more realistic case of the weakly anisotropic Heisenberg SG, by contrast, there no longer remains any 
global symmetry in fields. Hence, from symmetry, the situation is the same as that of the well-studied Ising SG. 
Meanwhile, in view of the fact that the Heisenberg SG possesses the nontrivial chiral degree of freedom which is 
totally absent in the Ising SG, the question of whether the ordering properties of the weakly anisotropic Heisenberg 
SG in fields are essentially the same as those of the Ising SG in fields seems not so trivial. This question is further 
promoted by the apparently contradicting experimental observations on the Ising-like and weakly anisotropic 
Heisenberg-like SGs^M 

In the present paper, we study both the spin-glass and the chiral-glass orderings of the weakly anisotropic 
Heisenberg SG in magnetic fields by extensive equilibrium Monte Carlo (MC) simulations^. We find a clear 
numerical evidence that a finite-temperature RSB transition occurs in the chiral sector, which also accompanies 
the simultaneous SG order. Thus, in spite of the similarity in the symmetry properties, the ordering properties of 
the weakly anisotropic Heisenberg SG model turn out to be quite different from those of the standard Ising SG. 
This highlights the importance of the chirality. 

The paper is organized as follows. In §2, we introduce our model and explain the details of the MC simulation. 
Various physical quantities calculated in the simulation are defined in §3, and the results of our numerical simulation 
are presented in §4. In §5, we perform the scaling analysis of the critical properties of the transition, and construct 
a phase diagram of the model in the temperature - magnetic field plane. Section 6 is devoted to summary and 
discussion. 



II. THE MODEL AND THE METHOD 

In this section, we introduce our model and explain some of the details of our numerical method. The model 
we consider is the isotropic classical Heisenberg model on a 3D simple cubic lattice defined by the Hamiltonian, 

N 

n = -Y / (j ij s i -s j + D^s ifi s jv )-Hj2s iz , (1) 

<ij> ^i,u—x,y,z i—1 

where Si = (Six, St y , Si Z ) is a three-component unit vector, and H is the intensity of magnetic field applied 
along the z direction. The isotropic nearest-neighbor exchange coupling Jy is assumed to take either the value 
J or — J with equal probability, while the nearest-neighbor random exchange anisotropy 's (/J,, v=x,y,z are 
spin-component indices) are assumed to be uniformly distributed in the range [-D : D], where D is the typical 
intensity of the anisotropy. We impose the relation = = D^j 1 . 

We perform equilibrium MC simulations on this model. In the present simulation, we fix D/J = 0.05, which 
is a typical value of D of real Heisenberg-like SG materials. Simulations are then performed for a variety of field 
intensities in the range H/J = 0.02 — 3.0. The lattices studied are simple-cubic lattices with N — L 3 sites with 
L = 4, 6, 8, 10, 12 and 16 with periodic boundary conditions. Sample average is taken over 64-800 independent 
bond realizations, depending on the system size L and the field intensity H. Limited amount of data are also 
taken for L = 20 in some cases (32 samples) to check the size dependence of physical quantities. 

To facilitate efficient thermalization, we combine the standard heat-bath method with the temperature-exchange 
technique^ The temperature-exchange trial is performed every heat-bath sweep. Typically, for the size L = 
16 (L — 20), we discard initial 8 x 10 5 (13 x 10 5 ) heat-bath sweeps and the temperature-exchange trials for 
equilibration, and use subsequent 8 x 10 5 (13 x 10 5 ) heat-bath sweeps and the temperature-exchange trials in 
calculating various physical quantities. Care is taken to be sure that the system is fully equilibrated. Equilibration 
is checked by the following procedures: First, we monitor the system to travel back and forth many times during the 
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TABLE I: Details of our MC simulations. H/ J represents the magnetic-field intensity, L the lattice size, N a the total 
number of samples, Nt the total number of temperature points used in the temperature-exchange run, T max /J and Tmm/J 
the maximum and minimum temperatures in the temperature-exchange run. 
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2~min / J 
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800 


26 


0.475 


0.085 
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800 


26 


0.475 


0.085 
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600 


26 


0.475 


0.113 
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10 


384 


42 
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0.115 




12 


256 


52 


0.40 


0.115 
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50 


0.35 


0.12 
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32 
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0.35 


0.1775 
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400 


26 


0.475 


0.085 




6 


400 


26 


0.475 


0.085 


0.5 
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400 


26 


0.475 


0.113 




10 


300 


42 


0.40 


0.115 




12 
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52 
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0.115 




16 
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0.12 
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0.1775 




4 


400 


26 


0.475 
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400 


26 


0.475 
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26 


0.475 


0.113 




10 


300 


42 


0.40 


0.115 




12 


256 


52 


0.40 


0.115 




16 


64 


50 


0.35 


0.125 



the temperature-exchange process (typically more than 10 times) between the maximum and minimum temperature 
points, and at the same time check that the relaxation due to the standard heat-bath updating is reasonably fast at 
the highest temperature, whose relaxation time is of order 10 2 Monte Carlo steps per spin (MCS). This guarantees 
that different parts of the phase space are sampled in each "cycle" of the temperature-exchange run. Second, 
we check the stability of the results against at least three times longer runs for a subset of samples. Error bars 
of physical quantities are estimated by the sample-to-sample statistical fluctuation over the bond realizations. 
Further details of our MC simulations are given in Table [I] 

III. PHYSICAL QUANTITIES 

In this section, we define various physical quantities calculated in our simulations. 

A. Chirality-related quantities 

We begin with the definition of the chirality. The local chirality at the z-th site and in the /i-th direction, Xi^i 
is defined for three neighboring Heisenberg spins by the scalar 

XiM = S i+ e. ■ (Si x Si.eJ, (2) 

where e M (/i = x, y, z) denotes a unit vector along the /i-th axis. By this definition, there are in total 3N local 
chiral variables in the system. The local chirality amplitude is then defined by 



^E E W . (3) 

i=l fi=x,y,z 

where (• • •} represents the thermal average and [• ■ •] represents the average over the bond disorder. The local 
chirality amplitude gives us the information of the extent of the non-coplanarity of local spin structures. 




By considering two independent systems ("replicas") described by the same Hamiltonian, one can define an" 
overlap of the chiral variable via the relation, 



1 N (a) (6) 

E E 

Xz/i X.i}_L 1 (4) 

i=l n=x,y,z 

where xl^ an d Xi b J represent the chiral variables of the replicas "a" and "b", respectively. In our simulations, we 
prepare the two replicas a and b by running two independent sequences of systems in parallel with different spin 
initial conditions and different sequences of random numbers. 

Since the present model does not possess any global symmetry, an odd quantity [(q x )] i s generally non-zero 
even in the hight-temperature phase. Taking this effect into consideration, the chiral-glass order parameter may 
be defined by 

?f } = [((Qx ~ [(Qx)]) 2 )} ■ (5) 
The associated chiral-glass susceptibility, normalized by the local amplitude x, is defined by 

A 

The Binder ratio of the chirality is defined by 
where 

Q [ x ] = [((Qx ~ [(Qx)]) 4 )} ■ (8) 

Here, g' x is normalized so that, in the thermodynamic limit, it vanishes in the high-temperature phase and gives 
unity in the ordered phase if the ordered state is non-degenerate. The distribution function of the chiral overlap 
q x is defined by 

Px(Q'x) = Ml'x-1x))} • (9) 

We consider the Fourier-transformed two-point chiral-glass correlation function C GG (fc) between the two local 
chiral variables in the /i-th and in the v-th directions. While one can define various types of correlation functions 
depending on the relative directions of the chiral variables (fi, v) and the direction of k, we consider here the 
parallel component C|l G (fe) where \i and v are both parallel with k. Here, we put k parallel with the x-direction, 
k = (k, 0, 0), so that n = v = x. Then, c|i G (fe) can be written in terms of the /c-dependent chiral overlap q x (k), 

4g(*) = [< \Qx(k)\ 2 >}, (10) 



1 N 

Qx( k ) = -fi E XixxfJ exp(ifc • ri), 



(11) 



where r, = (xi,yi,Zi) denotes the position vector of the chiral variable at the i-th site. The associated chiral 
correlation length, £ x , is defined by, 



2sin( 



c CG (o) 



(km/2) \ C CG (fe m ) 



- 1 



(12) 



where k m — (2w/L, 0, 0) is a wavevector of the minimum magnitude. 

In order to study the equilibrium dynamics of the model, we compute the equilibrium autocorrelation function 
of the chirality defined by 



<%(*) = (^E E X*(*o)Xi,.(* + *o)) 

. \ i=l fJ,=x,y,z I _ 



\(Qx)] 



(13) 



where the "time" t is measured in units of MCS. This chiral autocorrelation function is an "odd" quantity no"? 
invariant under the global flipping of the chirality. The corresponding "even" time correlation functions which is 
invariant under the global flipping of the chirality, may be defined by 



i—1 fj,=x,y,z 



(14) 



In displaying the data, we normalize these time correlation functions by their values at a unit time t = 1, i.e., we 
put 



c x (t) = c'(t)/cm, 



<rw/<T(i)- 



(15) 



Note that in the above definitions of the time-correlation functions (13) and (14), the second terms, [{q x }] and 
[(<?£)], have been subtracted, which are nonzero even in the high-temperature phase in the L — > oo limit due to 

~ (2) 

the absence of any global symmetry This subtraction guarantees that both C x (t) and q x (t) decay to zero as 

~ (2) 

t — > oo in the high-temperature phase. In the possible ordered phase, by contrast, both C x (t) and q x (t) decay 
to zero if the ordered state does not accompany the RSB, but tend to finite positive values if the ordered state 
accompanies the RSB. The latter property arises because, in the presence of RSB, the t — > oo limits of the first 
terms of Eqs.(13) and (14) are generally greater than the second terms, [{q x }] and [(</£)]• 

In computing the first terms of Eqs.(13) and (14), the simulation is performed according to the standard heat- 
bath updating without the temperature-exchange procedure, while the starting spin configuration at t = to is taken 
from the equilibrium spin configurations generated in our temperature-exchange MC runs. The second terms of 
Eqs.(13) and (14) are evaluated from the temperature-exchange MC runs. 



B. Spin-related quantities 

As in the case of the chirality, it is convenient to define an overlap variable for the Hciscnbcrg spin. In this case, 
the overlap might naturally be defined as a tensor variable q^ between the \i and v components (fx, v—x,y, z) of 
the Heisenberg spin, 



^ = ^f>l a) 4 6) . (m = w). ( 16 ) 



where and Sf^ are the i-th Heisenberg spins of the replicas a and b, respectively. 

In terms of these tensor overlaps, the "longitudinal" (parallel to the applied field) and the "transverse" (per- 
pendicular to the applied field) SG order parameters may be defined by 

5l 2) = K(«l - [(9l)]) 2 )], <?L=fe , (17) 
<?t } = [((?t - [(?t)]) 2 )], qr = ««* = + %v ■ (18) 

Note that as in the case of the chirality, the expectation value of the first moment has been subtracted. We also 
consider the spin-glass susceptibility for the transverse spin component defined by, 

XT = Nq { V (19) 

The longitudinal and the transverse Binder ratios are defined, respectively, by 
where 

i 4) - K(«l - [(a)]) 4 )] , (22) 

4 4) - [{(9r [<9T>]) 4 }] • (23) 



Here, g' h and g' T are normalized so that, in the thermodynamic limit, they vanish in the high-temperature phase 
and give unity in the ordered state if the ordered state is non-degenerate. 

The full spin-overlap distribution function may be defined in the tensor space with 3x3 = 9 components. Here, 
we consider the spin-overlap distribution function for the longitudinal and the transverse components, each defined 

by 

Wl) = [Wl-<7l)>] , (24) 

PM = [(6(4 r -< tr ))] , (25) 

where <7l and gx are defined by Eqs.(17) and (18), respectively. In the possible SG ordered state of the isotropic 
system, P s (<7t) develops a non-trivial shape in the thermodynamic limit due to the fact that the transverse- spin- 
overlap qt transforms non-trivially under the global spin rotation: See Refs. 38 and 39 for details. But here, the 
system is anisotropic so that no non-trivial structure arising from the uniform global spin-rotation is expected to 
arise in P s (gx) in the thermodynamic limit. 

We consider the Fourier-transformed spin-glass correlation functions both for the longitudinal and the transverse 
components, Cl(&) and Or(fe), which can be written in terms of the /c-dependent longitudinal and transverse spin- 
overlaps, qh(k) and qT(k), as, 



C L (fe) = [< |g L (fc)| 2 >], C T (k) = [< \q T (k)\ 2 >], (26) 



with 



N , N 



^ = 7j E ^M 6) ex P (*fc • n), 9t{k) = ^ £ S# • exp(zfc • n ). (27) 



i=l i=l 



where S^r = (Si X , Si y ) represents the transverse (xy) component of the Heisenberg spin. The associated longitu- 
dinal and the transverse spin correlation lengths, £l and £t, are defined by 



1 / ' , (28) 



2sin(fc m /2) V C L (fe m ) 



P 1 / Ct(0) 1 r9q , 

tT 2sin(fc m /2)Vc T (fc m ) ' 

respectively. 

The spin autocorrelation functions are defined both for the longitudinal and the transverse components by 

1 N 

c L (t) = ci(t)/ci(i), cut) = -]T[<MMMi + M}]-K9L)] , (30) 



i N 

C T (t) = C^(t)/<%(1), C^(*) = - ^[(Sff(to) • S<t(* + *o)>] - [(5t>] • (31) 



i=l 



The second terms, [(<7l)] and [(<Zt)], have been subtracted in the same context as in the definition of C x (t) for 
the chirality. These spin autocorrelation functions are computed in the same way as the chiral autocorrelation 
functions, and are normalized at their values at a unit time t = 1. The corresponding "even" time correlation 
functions can also be defined, though we skip their definitions here. 



IV. NUMERICAL RESULTS 



In this section, we show the results of our MC simulations on the anisotropic Heisenberg SG with D/J = 0.05. 



A. Time correlation functions of the chirality 
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In Fig2] we show the MC time dependence of the chiral autocorrelation function C x (t) for the size L = 16 
on log-log plots, for the fields (a) H/J = 0.05, (b) H/J = 0.5 and (c) H/J = 3.0, respectively. In the case 
Hj J = 0.05, we have checked that, in the time range shown, the data can be regarded as those of the bulk, since 
no appreciable size effect is discernible between the data of L = 16 and L = 20, the latter data being shown with 
lines in the figure. In Fig. 2, we show the MC time dependence of the corresponding even quantities, q x \t), for the 
fields (a) H/J = 0.05, (b) H/J = 0.5 and (c) H/J = 3.0, respectively. As can clearly be seen from these figures, 

~ (2) 

in the investigated time range, the observed behavior of C x (t) is essentially the same as that of q x '(t). 

For the field H/J = 0.05, as can be seen from the Figs. 1(a) and 2(a), C x (t) and q x \t) exhibit either a down- 
bending or an up-bending behavior depending on whether the temperature is higher or lower than a borderline 
value T/J ~ 0.21. while just at this borderline temperature a straight-line behavior corresponding to a power-law 
decay is observed. This indicates that the chirality exhibits a phase transition into the low-temperature ordered 
phase where the replica symmetry is spontaneously broken. In the case of H/J — 0.5, as can be seen from Figs. 1(b) 
and 2(b), the data at T/J 5 0.24 show a slight up-bending tendency at a short time t ~ 10 2 while the data at 
T/J c 0.21 show a gradual down-bending tendency at longer times t ~ 10 3 , and the chiral transition temperature 
appears to lie somewhere between T/J ~ 0.21 and 0.24. The transverse spin autocorrelation, which is to be shown 
in Fig. 3(c) below, however, suggests the transition temperature T g ~ 0.20 — 0.21. Then, we finally estimate the 
transition temperature T g /J = 0.21 ± 0.02 for H/J = 0.5 (see below). In the case H/J = 3.0, by contrast, C x (t) 

(2) 

and q x (t) always exhibits a down-bending behavior in the temperature range studied, suggesting the absence of 
a phase transition, at least in the temperature range T/J > 0.125. 

~ ( 2) 

Since very much similar behaviors are observed in C x (t) and in q x (t) , we shall show in the following subsections 
the data of the autocorrelation function C (t) only. 



B. Time correlation functions of the spin 

In FigsOl we show for the fields H / J — 0.05 and 0.5 the MC time dependence of the autocorrelation functions of 
the transverse component of the spin Cr(i), Figs, (a) and (c), and those of the longitudinal component of the spin 
Cl, Figs.(b) and (d). Finite-size effect evaluated from the difference between the L — 16 and L — 20 data turns out 
to be rather small (but not completely negligible) in the longitudinal component, whereas it is more appreciable 
in the transverse component. Although a finite-size effect is not completely negligible here, the autocorrelation 
functions of the spin for H/J = 0.05 turn out to behave quite similarly to those of the chirality. Namely, both 
Ch{t) and Ct(£) exhibit either a down-bending or an up-bending behavior depending on whether the temperature 
is higher or lower than the borderline value T/J ~ 0.21, while just at this borderline temperature a straight-line 
behavior corresponding a power-law decay is observed. This indicates that the spin exhibits an RSB transition at 
the same temperature where the chirality exhibits an RSB transition. The simultaneous occurrence of the spin 
and the chirality orderings is quite natural in the presence of the random anisotropy. 

For the field H/J — 0.5, the transverse component Ct(£) shown in Fig. 3(c) exhibits a rather clear up- 
bending/down-bending behavior with the borderline temperature T/J ~ 0.20 — 0.21, which might be com- 
pared with the less clear behavior of C x [t) where the transition temperature appears to lie somewhere between 
T/J ~ 0.21 — 0.24. Since the transition here is expected to be a simultaneous spin and chiral transition, we now 
estimate the transition temperature of the field H/J — 0.5 to be T g = 0.21 ±0.02. Meanwhile, the longitudinal au- 
tocorrelation function C^(t) for H/J = 0.5 exhibits a much noisier behavior as shown in Fig. 3(d). This is because, 
in large fields, the second term of Eq.(30), [(<7l)], becomes large, while the longitudinal spin correlation function 
itself becomes small in magnitude which is obtained as a difference between the two large numbers. Note that the 
up-bending behavior is discernible at a short time t ~ 10 2 , which is also observed in the chiral autocorrelation in 
a less pronounced manner. In any case, due to the noisiness, it seems not possible to identify the transition point 
from Fig. 3(d). 

From the behaviors of the chiral and the spin autocorrelation functions shown above, we conclude that the spin 
and the chirality exhibit an RSB transition simultaneously at a finite temperature T g j J = 0.21 ± 0.02 for both 
the fields H/J = 0.05 and 0.5. The estimated transition temperature T g /J = 0.21(2) is lower than the transition 
temperature of the corresponding zero-field model with the same magnitude of anisotropy, i.e., T g /J ~ 0.24j 2 ^* 2 ^ 
indicating that applied magnetic fields suppress the spin-glass (chiral-glass) ordering for weaker fields. Indeed, this 
is exactly the feature expected for the experimental AT-linc. Interestingly, the estimated transition temperature 
under fields, T g / J ~ 0.21, comes very close to the chiral-glass transition temperature of the fully isotropic model 
in zero field, which was estimated to be Tqq/ J ~ 0.20i22i This is exactly the feature expected for the experimental 
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FIG. 1: (Color online) Temporal decay of the autocorrelation function of the chirality C x (t) defined by Eqs. (13) and (15), 
for the fields (a) H/J = 0.05, (b) H/J = 0.5 and (c) H/J = 3.0. The lattice size is L = 16. The data at T = T g are 
given in red (by filled symbols). In Fig. (a), in order to check the finite-size effect, the data of L = 20 are plotted with lines 
at temperatures T/J — 0.18-0.24 with an interval of 0.01. The estimated transition temperature is T g /J ~ 0.21 both for 
H/J = 0.05 and H/J = 0.5. 



GT-linc. Hence, the suppression of T g (H) due to weaker fields (AT-line) as well as the robustness of it with 
respect to stronger fields (GT-line) are consistent with the experimental observation for the weakly anisotropic 
Heisenberg-like SGsi^ 



C. Comparison with the autocorrelation functions of the 3D Ising SG in fields 

For comparison, we also calculate the autocorrelation function of the spin Cl (t) defined by Eq. (|30[) for the 3D 
Ising SG with the ± J coupling for the field H/ J = 0.05. Note that the Ising SG in fields shares the same symmetry 
property as the anisotropic Heisenberg SG in fields, i.e., the absence of any global symmetry. Nevertheless, the 
Ising SG does not possess any chiral degree of freedom. Thus, the question of whether the 3D Ising SG behaves 
either similarly or differently from the weakly anisotropic Heisenberg SG would be of special interest. 

The data of Cl(£) of the Ising SG are shown in Fig0|for the size L — 20. As can be seen from figure, the behavior 
of Cjjit) of the Ising SG differs significantly from that of the weakly anisotropic Heisenberg SG. Although the 
temperature range studied is as low as about 60% of the zero- field transition temperature T g {H = 0) / J ~ 1 .1 , which 
is expected to be deep in the ordered state according to a tentative estimate of Ref. 18 , no clear up-bending behavior 
as observed in the weakly anisotropic Heisenberg SG is observed here. Instead, (t) and Cl(£) persistently exhibit 
an almost linear behavior even at the lowest temperature studied T/J = 0.6. This is illustrated in Fig. 4(b) where 
the Ci,(t) data at the two lowest temperatures studied are shown. A comparison of the L — 20 data with the L = 16 
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FIG. 2: (Color online) Temporal decay of the time-correlation function of the chirality <jx\t) defined by Eqs.(14) and (15), 
for the fields (a) H/J = 0.05, (b) H/J = 0.5 and (c) H/J = 3.0. The lattice size is L = 16. The data at T = T g are 
given in red (by filled symbols). In Fig. (a), in order to check the finite-size effect, the data of L = 20 are plotted with lines 
at temperatures T/J — 0.18-0.24 with an interval of 0.01. The estimated transition temperature is T g /J ~ 0.21 both for 
H/J = 0.05 and H/J = 0.5. 



data indicates that some amount of finite-size effect still remains. Nevertheless, an almost linear behavior without 
any discernible up-bending tendency is robustly observed in common both for L = 16 and L = 20, suggesting that 
this feature is a bulk property. 



D. The Binder ratio of the chirality 



In FigEl we show the temperature- and size- dependence of the chiral Binder ratio g' x for the field H/J = 0.05. 
The data of g' x has a negative dip at a size-dependent temperature T = Tdi p (L), which, with increasing L, tends to 
deepen and shift to lower temperatures. The existence of a persistent negative dip with increasing depth, is a sign 
of a phase transition occurring at T = T<np(L = oo). By extrapolating Tdip(^) to L — oo, as shown in the inset 
of FigO the bulk chiral-glass transition temperature is estimated to be T g /J = 0.22(2). The estimated transition 
temperature agrees well with the corresponding estimate based on the autocorrelation functions given above. 

As argued in RefiSL^LIS for the isotropic case, the existence of a persistent negative dip is a sign of a one-step- 
like RSB transition. In the present anisotropic model under fields, there no longer exists a global spin-reflection 
symmetry in contrast to the isotropic model under fields. Hence, the behavior expected for g' here might be 
the one for the one-step RSB system without a reflection symmetry. Such a system was theoretically analyzed in 
Ref4i, where the behavior of the Binder ratio g' x in the thermodynamic limit was reported as shown in FigHjJ 41 
Indeed, the overall behavior of our present g' x shown in Fig[B] seems consistent with such a behavior. 
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FIG. 3: (Color online) Temporal decay of the spin autocorrelation functions: (a) the transverse one Cx(i) for H/J = 0.05, 
(b) the longitudinal one C^it) for H/J = 0.05, (c) the transverse one for H/J = 0.5, and (d) the longitudinal one for 
H/J = 0.5. The lattice size is L — 16. The data at T = T a are given in red (by filled symbols). In the figures, in order to 
check the finite size effect, the data of L = 20 are plotted with lines at temperatures T/J — 0.18-0.24 with an interval of 
0.01. The estimated transition temperatures is T g /J ~ 0.21 both for H/J — 0.05 and H/J — 0.5. 



E. The Binder ratio of the spin 

In FigUfa), we show the temperature and size dependence of the Binder ratio of the longitudinal component 
of the spin g' h for the field H/J = 0.05. Although g' h for smaller sizes L < 12 does not show any characteristic 
feature suggestive of a phase transition, the one for larger sizes L > 16 tends to exhibit a negative dip similar 
to the one as observed in the chiral Binder ratio g' . Then, for large enough L, the overall shape of g'^ would be 
similar to that of g' x , although the one for smaller L is very different. 

In Fig. 7(b), we show the temperature and size dependence of the Binder ratio of the transverse component of 
the spin g' T for the field H/J = 0.05. Although g' T for smaller sizes L < 12 exhibits a single maximum with a 
crossing point occurring around T/J ~ 0.28, the one for larger sizes L > 16 exhibits double maxima. The apparent 
crossing point observed for smaller sizes at T/J ~ 0.28 disappears for larger lattices, indicating that it does not 
correspond to a true phase transition point. Among the two maxima of g' T observed for L > 16, the one at a 
higher temperature is gradually suppressed with increasing L, while the one at a lower temperature tends to be 
enhanced. Presumably, in the L —> oo limit, the peak at a higher temperature will disappear, while the peak at a 
lower temperature will survive. Then, for large enough L, the overall shape of g' T would be similar to that of g' 
and <7l, although the one for smaller L is very different. Our observation here that the Binder ratio of the chirality, 
g' x , and those of the spin, g' L and g' Tl asymptotically show mutually similar behavior, which resembles the one 
depicted in Fig. 6, seems consistent with the spin-chirality decoupling-recoupling scenario. Indeed, the recoupling 
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FIG. 4: Temporal decay ol the spin autocorrelation function, denned by Eq. IMOt . of the ±J 3D Ising SG in a magnetic 
field of 77/ J = 0.05. The system size is L — 20 averaged over 100 samples. In (a), the data are plotted in the temperature 
range T/J — 0.60 — 1.3 with an interval of 0.1. The data of the two lowest temperatures T/J — 0.60 and 0.70 are shown 
in (b), together with the fitted straight lines. Even at the lowest temperature studied, no up-bending behavior is observed. 
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FIG. 5: (Color online) The temperature and size dependence of the Binder ratio of the chirality g' x for a field 77/ J = 0.05. 
In the inset, we show the reciprocal size dependence of the dip temperature of g' x . Dashed line represents a linear fit of the 
data. The bulk chiral-glass transition temperature is estimated to be T g /J = 0.22(2). 



length-scale estimated in Refi 2 V 8 , L x ~ 20, is consistent with the observed behavior. 



F. Overlap distribution function of the chirality 



We show in Figs |Hl the size dependence of the overlap distribution function of the chirality, P x (q x ), for several 
cases, i.e., (a) H/J = 0.05 and T/J = 0.18, (b) H/J = 0.5 and T/J = 0.19, (c) H/J = 0.5 and T/J = 0.12, and 
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FIG. 6: A sketch of the typical temperature dependence of the Binder ratio g' in the thermodynamic limit, which is 
expected in a system without a reflection symmetry exhibiting a one-step RSB transition. 
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FIG. 7: (Color online) The temperature and size dependence of the Binder ratio of (a) the longitudinal component of the 
spin g' L , and of (b) the transverse component of the spin g' T , for a field H/J — 0.05. 



(d) H/J = 3.0 and T/J = 0.13. 

In the case H/ J = 0.05, in addition to a primary peak corresponding to q x — q^ A > 0, which grows and sharpens 
with increasing L, there appears a second peak at around q x = 0, which also grows and sharpens with increasing 
L. The existence of two distinct peaks, both growing and sharpening with increasing L, is a clear indication of 
the occurrence of RSB. As reported in Refi21, P x (q x ) in zero field exhibits a feature of a one-step-like RSB, i.e., 
a central peak at q x — coexisting with self-overlap peaks at q x — ±q x A - The P x {q x ) observed here may be 
regarded as the in- field counterpart of the zero-field P x (q x ) with a feature of such a one-step-like RSB. Indeed, if 
one closely looks at P x {q x ) shown in FigOJa), one sees that a broad peak, which is a remnant of the q x = —q x A 
peak of the zero-field model, is discernible for smaller sizes L < 8, whereas, for larger sizes, this q x = —q x A peak 
disappears and the q x ~ peak begins to grow. Interestingly, for the size L = 8, P x {q x ) possesses three broad 
peaks, at around q x — ±q x A and q x = 0. Such a three-peak structure is rarely seen in a system exhibiting the full 
RSB, and gives a further indication that the RSB occurring here is the in- field counterpart of the one-step-likc 
RSB. 

For H/J — 0.5, P x (q x ) shows a similar behavior as that for H/J = 0.05. Namely, for larger L, it exhibits 
two distinct peaks, both growing and sharpening with increasing L, whereas it exhibits three broad peaks for an 



intermediate L. As compared with the H/J = 0.05 case, the second peak is located slightly off q x = 0, reflectiAg 1 
the fact that the higher field breaks the q x <-> — q x symmetry more strongly. Anyway, our data shown in Figs[SJa)- 
(c) give a strong numerical support that there indeed occurs a chiral-glass transition at a finite temperature and 
that the chiral-glass ordered state accompanies a one-step-like RSB. 

For the still higher field H/J = 3.0, in contrast to the cases of H/J = 0.05 and 0.5, the double-peak behavior 
of P x (q x ) is not observed even at the lowest temperature studied T/J = 0.13. Thus, for H/J = 3.0, no sign of 
RSB transition is observed down to this low temperature. 




FIG. 8: (Color online) The chiral-overlap distribution functions for the field and the temperature, (a) H/J = 0.05 and 
T/J = 0.18, (b) H/J = 0.5 and T/J = 0.19, (c) H/J = 0.5 and T/J = 0.12, and (d) H/J = 3.0 and T/J = 0.13. The 
transition temperature is T g / J ~ 0.21 for both cases of H/J = 0.05 and 0.5, while it is lower than 0.13 for H/J = 3. 



G. Overlap distribution function of the spin 



We show in FigsEjlthe size dependence of the overlap distribution function of (a) the longitudinal component 
of the spin P s (q^), and of (b) the transverse component of the spin P s (qT), for the field H/J — 0.05 and at a 
temperature T/J = 0.18. The longitudinal spin-overlap distribution function P s (qh) for smaller sizes L < 16 
possesses only a single growing peak at around <7l — 0.2, in apparent contrast to the double-peak structure 
observed in the chiral-overlap distribution. Quite interestingly, however, for the largest size studied L = 20, one 
sees that the second peak just begins to emerge at around (?l — 0.1 (see the arrow in Fig. 9(a)). This second 
peak is reminiscent to the one observed in the chiral-overlap distribution function of Fig. 8, though it appears 
here only for the largest size in a less pronounced manner: It appears to be an echo of the strongly diverging 
q x ~ peak observed in the chiral-overlap distribution. Our observation that the one-step RSB-like structure of 



the overlap distribution appears in the chiral sector from smaller sizes, while it appears in the spin sector in a 
pronounced manner only for larger sizes, suggests that the order parameter of the present one-step- RSB transition 
might be the chirality, rather than the spin. Again, this observation is fully consistent with the spin-chirality 
decoupling-recoupling scenario i 24 ' 2 ^ 

As shown in Fig. 9(b), the behavior of the transverse-spin-overlap distribution function P s (qT) is more complex. 
Such a behavior is certainly consistent with the occurrence of an RSB transition. 
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FIG. 9: (Color online) The overlap distribution functions of 
transverse component of the spin, for a field H/ J = 0.05 at 
Tg/J ~ 0.21. In Fig. (a), a sign of the second peak just begins 
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(a) the longitudinal component of the spin, and of (b) the 
a temperature T/J — 0.18. The transition temperature is 
to emerge for the largest size L — 20 at the arrow position. 



H. Correlation lengths 

We show in Figs^|the temperature and size dependence of the normalized correlation lengths of (a) the chirality 
£ X /L, of (b) the longitudinal component of the spin £l/L, and of (c) the transverse component of the spin £t/L, 
for the held H/J — 0.05. The location of T g (H) obtained from the autocorrelation functions are displayed with 
arrows in the figures. 

Concerning the chirality, as can be seen from the inset of Fig. 10(a). a clear crossing behavior is observed at 
T/J ~ 0.22 for larger sizes L > 12. The crossing temperature T/J ~ 0.21 observed for larger sizes turns out to 
be close to our previous estimate of T g . 

Concerning the longitudinal and the transverse components of the spin, a crossing is observed at T g J J ~ 0.21 for 
smaller sizes L < 8. For L > 10, however, £/L decreases with L at any temperature studied, no longer exhibiting 
a crossing at T g / J ~ 0.21. For even larger sizes L = 16 and L = 20, a tendency of crossing re-appears at around 
T g /J ~ 0.21, but now at a lower value of £/L. 

Such a complex size dependence of the spin correlation lengths may naturally be interpreted by the spin-chirality 
decoupling-recoupling scenario in the following way: For smaller sizes, i.e., at shorter length scales, the spin and 
the chirality are trivially coupled in correlations irrespective of the anisotropy, so that the crossing in the sizes 
L < 8 may reflect this trivial coupling at short length scales. For larger sizes, i.e., at long length scales, the spin is 
"recoupled" to the chirality via the anisotropy in a way different from the trivial coupling at short length scales. 
Thus, the two different types of crossing is expected in the normalized correlation lengths for smaller and for larger 
sizes, which is exactly the behavior observed in Figs. 10(b) and (c). The characteristic length scale separating the 
coupling and the recoupling regimes was estimated to be about 20 lattice spacingsi 29 ' 38 Therefore, although we 
cannot simulate here the sizes larger than L = 20 due to the lack of our computational capability, we do expect 
that a clear crossing behavior will eventually set in for L > 20. 



V. THE CRITICAL PROPERTIES AND THE PHASE DIAGRAM 



In this section, we further analyze the nature of the chiral-glass (spin-glass) transition observed in the previous 
section by means of a scaling analysis, and construct a magnetic phase diagram of the model in the temperature 
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FIG. 10: (Color online) The temperature and size dependence of the normalized correlation lengths of (a) the chirality, of 
(b) the longitudinal component of the spin, and of (c) the transverse component of the spin, for a field H/J = 0.05. In 
Fig. (a), a magnified figure is given in the inset for the three largest sizes, L = 12, 16 and 20. The arrows in the figures 
indicate the locations of T g (H), i.e., T g /J = 0.21. 



- magnetic field plane. 

We apply a dynamical scaling analysis to the chiral autocorrelation functions in order to estimate the critical 
exponents of the transition. The standard bulk dynamical scaling form is assumed for the autocorrelation functions, 
C x (t), C L (t) and C T (t), 

C(t) « |(T- T g )/jff(t\(T - T g )/JD (32) 

where T g is the transition temperature determined in the previous section, while (3, v and z refer to the order 
parameter, the correlation-length and the dynamical exponents, respectively. The exponents f3 and zv are to be 
determined so that a good data collapse is obtained in the scaling plot. The quality of the scaling plot is judged 
by eyes. 



A. Dynamical scaling analysis of the chiral autocorrelation function 

In Fies llll we show the scaling plots of the chiral autocorrelation function C x (t) for the fields (a) H/J = 0.05 
and (b) H/J = 0.5. The transition temperature is fixed to be T CG /J = 0.21 both for H/J = 0.05 and H/J = 0.5, 
as determined in the previous section. 

As can be seen from the figures, the chiral autocorrelation function scales well both above and below T g , with 
Pcg = 0.8 and zqgvqg = 5.0 for H/J = 0.05, and with /3cg = 0.7 and zcg^cg = 5.5 for H/J — 0.5. Examining 
by eyes the quality of the scaling plots with varying the assumed exponents values, we finally quote /3cg = 0.8 ±0.2 



and zcg^cg = 5.0±1.0 for H/J = 0.05, and /3 CG = 0.7±0.2 and z CG ^ CG = 5.5±1.0 for H/J = 0.5. The estimated 
chiral-glass exponents are not far from the corresponding zero-field values /?cg — 1-1— and zcg^cg — 4. 5^ 
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FIG. 11: (Color online) Dynamical scaling plots of the chiral autocorrelation function C x (t), for fields (a) H/J — 0.05 and 
(b) H/J = 0.5. 



B. Dynamical scaling analysis of the spin autocorrelation functions 

Next, we apply a dynamical scaling analysis to the spin autocorrelation functions. In Figs^J we show the scaling 
plots both for (a) the longitudinal and for (b) the transverse components of the spin for the field H/J = 0.05. 
The transition temperature is fixed to be T g /J = 0.21. In Ref. 3 -, we observed that the quality of the scaling plot 

of the spin time-correlation function was not as good as its chiral counterpart q% ■ . As argued in Ref^, if 
the chirality, rather than the spin, is the order parameter of the transition, in the time regime shorter than the 
recoupling time scale (estimated to be 10 5 ~ 10 6 MCS), the spin time-correlation might not reach an asymptotic 
scaling regime even when the chiral time-correlation reaches an asymptotic scaling regime. Hence, although the 
scaling of Cx(i) and Cr(£) turns out to be fairly good as shown in Fig. 12, caution might be required here in 
regarding the fitted values of /?sg and .zsg^SG as true asymptotic spin exponents, since our observation time 
window is shorter than the expected recoupling time scale. 



C. Chiral-glass and spin-glass susceptibilities 

In order to estimate the chiral-glass and the spin-glass susceptibility exponents 7cg and 7SG ; we plot in Figs ll3l 
(a) the reduced chiral-glass susceptibility x x , and (b) the transverse spin-glass susceptibility xt, as a function of 
the reduced temperature (T — T g )/T g for the field H/J — 0.05 on a log-log plot. The chiral-glass and the spin-glass 
susceptibilities in the thermodynamic limit are expected to behave as 

X X ~\T-T g \'^, XT~\T~T g \-^. (33) 

The T g value is fixed to be T g /J = 0.21. From the asymptotic slope of the data, the exponents 7cg and 7sg are 
estimated to be 7cg = 2.1(2) and 7sg = 2.4(2), respectively. Here, the obtained values of 7sg and 7cg turn out 
to be rather close. Indeed, the spin-chirality decoupling-recoupling scenario predicts an equality 7sg = 7cg- 



D. H-T phase diagram 

Finally, we show in Fig^]a phase diagram of the present model with D/J = 0.05 in the temperature-magnetic 
field plane. The obtained phase diagram has two noticeable features. First, for weaker fields H/J < 0.05 = D/J, 
the transition line T g (H) exhibits an AT- like like behavior, i.e., the spin-glass (chiral-glass) transition temperature 
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FIG. 12: (Color online) Dynamical scaling plots of the autocorrelation functions of (a) the longitudinal component of the 
spin Cl, and of (b) the transverse component of the spin Cl, for a field H/ J = 0.05. 




FIG. 13: (Color online) Temperature and size dependence of (a) the chiral-glass susceptibility, and of (b) the transverse 
spin-glass susceptibility, for a field H/J = 0.05. The transition temperatures is set to be T g /J = 0.21. 



is suppressed rapidly under weak applied fields. In sharp contrast to this, for stronger fields of H/ J ;J 0.05 = D/J, 
the spin-glass (chiral-glass) transition temperature turns out to be rather insensitive to the field intensity, indicating 
that the spin-glass (chiral-glass) ordered state is robust against applied fields. The estimated T g (H) value stays 
almost unchanged between the field values H/J = 0.05 and 0.5, the latter value is already more than twice 
the zero-field transition temperature. Indeed, the T g (H) value in this field range is very close to the zero-field 
chiral-glass transition temperature of the fully isotropic model. This is exactly the feature of the GT line. 

Overall, the obtained phase diagram is consistent with the experimental phase diagram of Heisenberg-like SG 
magnetsi^ In order to make a more quantitative comparison, however, we have to determine the phase boundary 
in the lower field regime H/J < 0.05 more precisely, and to map out a phase diagram with varying the anisotropy 
strength D. This is a computationally demanding task, and is left as a future task. 
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FIG. 14: (Color online) The phase diagram of the anisotropic Heisenberg SG with the random anisotropy of D/J — 0.05 in 
the temperature- magnetic field plane. "CG" , "SG" and "PM" stand for the "chiral-glass" , "spin-glass" and "paramagnetic" 
phase, respectively. 



VI. SUMMARY AND DISCUSSION 



In summary, we performed a large-scale equilibrium Monte Carlo simulation of the weakly anisotropic 3D 
Heisenberg SG in magnetic fields, paying attention both to the spin and to the chirality. The model is expected 
to be a reasonably realistic model of many real Heisenberg-like SG magnets. Due to the presence of both the 
random anisotropy and the magnetic field, the model lacks in any global symmetry so that the transition, if any, 
is expected to be a pure RSB one. Among other things, we have found a clear numerical evidence that such a pure 
RSB transition indeed occurs at a finite temperature simultaneously in the spin and in the chiral sectors. The 
temporal decay of the time correlation functions of the chirality and of the spin, as well as the overlap distribution 
function of the chirality and the spin, gave a strong numerical evidence of the occurrence of a thermodynamic RSB 
transition in applied fields. We feel this finding for the weakly anisotropic Heisenberg SG to be rather remarkable, 
in view of the long controversy in the community concerning the in-field ordering properties of the Ising SG, which 
apparently looks much simpler. 

Indeed, a comparison with the 3D Ising SG in fields have revealed that the weakly anisotropic Heisenberg SG 
in fields behaves differently from the Ising SG in fields, a sign of the RSB transition being much clearer in the 
weakly anisotropic Heisenberg SG than in the Ising SG. The observed stronger sign of the RSB transition is closely 
connected to the one-step-like nature of its RSB pattern, which is certainly not the case in the Ising SG. Such a 
difference in the ordering properties of the two models might at first sound surprising, since the Ising SG in fields 
and the weakly anisotropic Heisenberg SG in fields share the same symmetry property. However, in spite of the 
similarity in symmetry, a significant difference exists in the two models in that the Ising SG does not possess any 
chirality degree of freedom in contrast to the weakly anisotropic Heisenberg SG. Then, our observation highlights 
the possible importance of the chirality in realizing the RSB transition in the weakly anisotropic Heisenberg SG. 
This gives a strong support to the chirality scenario of the SG transition of Heisenberg-like SGsi^iS^ 

We also have found that the sign of the RSB transition is often observed more strongly in the chirality-related 
quantities than in the spin-related quantities. For example, the second peak of the overlap distribution function 
shows up clearly even from smaller lattices for the chirality, but more weakly and only for larger lattices for the 
spin: Compare Fig. 8(a) and Fig. 9(a). This observation seems consistent with the chirality scenario that the order 
parameter of the transition is the chirality whereas the spin is recoupled to the chirality beyond a certain crossover 
(recoupling) length scale of about 20 lattice spacings. We have also observed that the behaviors of certain spin- 
related quantities, e.g., the Binder ratio or the normalized correlation length, differ between in smaller and in 
larger lattices, the borderline size being around L ~ 20. Again, these behaviors could naturally be explicable on 



the basis of the spin-chirality decoupling-recoupling scenario i^^S^ ^ 

We also have constructed a rough H-T phase diagram of the model. It has turned out that T g (H) is rapidly 
suppressed by weak applied fields, H IS D, but is then kept almost unchanged up to rather higher fields of 
H/J = 0.5 1=3 2T g (0). The phase boundary in the weak field regime resembles the AT-line, while that in the 
high field regime resembles the GT line, although its nature is very different from the mean-field AT- or GT- 
line i 5 i 6 i 7 These features of the phase diagram are consistent with the experimental result of the weakly anisotropic 
Heisenberg SG<^ The obtained phase diagram is also fully consistent with the one expected from the chirality 
scenario Indeed, the chirality scenario predicts that, in the low-field regime where the anisotropy is 
important, the change in the broken symmetry from the zero-field case causes a crossover, accompanied with a 
rapid suppression of the transition temperature, while in the high-field regime where the anisotropy is negligible 
relative to the applied magnetic field, the SG transition line should essentially be given by the chiral-glass transition 
line of the fully isotropic system^ 

Another remarkable feature of the present phase diagram is that the GT-like phase boundary is quite robust 
against magnetic fields, extending without much reduction to the field at least as twice large as the zero-field 
transition temperature. Such a robustness of the SG ordered phase under fields is also in accord with the chirality 
scenario, since the chirality, regarded as a hidden order parameter of the SG transition, is only weakly coupled 
to the magnetic field. Note that the magnetic field couples directly to the spin via the Zeeman term in the 
Hamiltonian, not to the chirality, while the effective coupling between the field and the chirality is rather weak. 

Thus, overall, our observations on the weakly anisotropic Heisenberg SG gives a strong support to the spin- 
chirality decoupling-recoupling scenario In other words, it would be difficult to interpret all of our present 
observations based on the more conventional view that the spin is the order parameter of the transition and the 
chirality is just a composite (secondary). 

A decisive test of the spin-chirality decoupling-recouping scenario would be to directly probe the chiral-glass 
order. Recently, a proposal was made concerning the way to measure the linear and the nonlinear chiral suscep- 
tibilities of canonical SGs by using the Hall probe, together with the scaling prediction from the spin-chirality 
decoupling-recoupling scenario42^ Interestingly, very recent experiments have given support to these predictions 
from the chirality scenario ^MJ^ 

In concluding this paper, we wish to emphasize again that, even if one sets aside the question of the de- 
tailed mechanism of the SG transition or the validity of the chirality scenario, the present numerical results have 
demonstrated beyond reasonable doubt that the weakly anisotropic Heisenberg SG exhibits a thermodynamic RSB 
transition in applied magnetic fields. In view of the fact that the property of the apparently simpler, strongly 
anisotropic Ising SG still remains unclear in spite of the long controversy, this finding seems surprising and some- 
what ironical. Clarifying the relation between the weakly anisotropic Heisenberg SG possessing the chiral degree 
of freedom and the strongly anisotropic Ising SG not possessing the chiral degree of freedom is still open, and is 
left for future studies. 
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